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Abstract 

We write a first order symmetric hyperbolic system coupling the 
Riemann tensor with the dynamical acceleration of a prefect relativis- 
tic fluid. We determine the associated, coupled, Bel - Robinson type 
energy, and the integral equality that it satisfies. 

Systeme de Bianchi-Euler pour un fiuide relativiste, et energie 
de type Bel-Robinson. 

Resume. On ecrit un systeme symetrique hyperbolique satisfait par le 
tenseur de Riemann de l'espace temps et l'acceleration dynamique d'un fiuide 
parfait relativiste. On determine l'energie du type Bel-Robinson correspon- 
dante, et l'egalite integrate qu'elle satisfait. 



1 Introduction. 



The effective strength of the gravitational field lies in the Riemann tensor of 
the spacetime metric. Its evolution is governed by the so - called higher order 
equations (Bel 1958, Lichnerowicz 1964), deduced from the Bianchi identi- 
ties. The system satisfied by the trace free part of the Riemann tensor, the 
Weyl tensor, was some time ago recognized as a linear, first order symmetric 
hyperbolic system (FOSH), with constraints, homogeneous in vacuum. See 
H. Friedrich (1996) and references therein. The evolution equations for the 
Riemann tensor itself have also been written (CB-Yo 1997) as a FOSH sys- 
tem, made explicit in terms of four two- tensors, introduced by Bel (1958) 
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in a general setting, the electric and magnetic gravitational fields and cor- 
responding inductions relative to a Cauchy adapted frame. In the presence 
of sources the higher order equations, now renamed Bianchi equations, are 
no longer homogeneous; their right hand sides are linear in the covariant 
derivative of the stress energy tensor of the sources. 

In this article we consider the case of perfect fluid sources. The Einstein 
equations with such sources have long ago (CB 1958) been proved to be a well 
posed Leray - hyperbolic system (with constraints). The fluid equations have 
also been written as a FOSH system (in special relativity, K.O. Friedrichs 
1969; in general relativity Ruggeri and Strumia 1981, Anile 1982, and Rendall 
1992 who proved a theorem valid also for isolated fluid bodies with special 
equations of state). It seemed however interesting to have a system of equa- 
tions which would be a FOSH system both for the gravitational field, namely 
the Riemann tensor of space time, and the fluid variables. Such a system has 
been written in lagrangian variables (that is, in a frame whose timelike axis 
is tangent to the fluid flow lines) by H. Friedrich 1998, who used the Weyl 
tensor and by CB-Yo 2001, using directly the Riemann tensor. In this paper 
we use eulerian variables, that is a frame adapted to the usual 3+1 slicing 
of the space time, with time axis orthogonal to the space slices, which we 
call a Cauchy adapted frame. We obtain a FOSH system for the Riemann 
tensor and the dynamical fluid acceleration. The energy corresponding to 
the FOSH system that we obtain, and call the Bianchi - Euler system, is 
the sum of the usual Bel Robinson energy of the gravitational field, and the 
dynamical acceleration energy of the fluid. It is not conserved in general - 
no more than the Bel - Robinson energy in vacuum - but its evolution can 
be controlled through the equations we obtain. 

2 Einstein equations with fluid sources. 

The Einstein equations with source a stress energy tensor are: 

Ra/3 = Pa/3, Pa/3 = T a/3 — -QapT 

The stress energy tensor of a perfect relativistic fluid is 

T a(3 = (p + p)u a u p + pg a/3 

with u a the unit kinematical velocity satisfying u a u a — — 1; and p,p,S the 
specific energy, pressure and entropy. These thermodynamic quantities are 
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assumed to be positive. They are linked by an equation of state: 

H = IM(p,S) (2.1) 
The enthalpy index / of the fluid is given by the identity 

f(p,S)= X ew / P 7 %, . (2-2) 
fj,(p,S)+p 



where x is a constant such that / 2 has the dimensions of an energy density^. 

The dynamical velocity C a , a vector tangent to the flow lines, incor- 
porates information on the kinematic velocity u a and the thermodynamic 
quantities. It is defined by 

C a = fu a , (2.3) 

hence it satisfies the relation 

C a C a = -f. (2.4) 

In terms of the dynamical velocity the equations of a perfect fluid in an 
arbitrary frame are 

C a V a C p + fdj = C a (V a C p - VpC a ) = (2.5) 

and, with fj,' p (p, S) = dp,(p, S)/dp, a given function of p and S, 

V a C a + (p> p - l)_-V a C^ = (2.6) 



C *y a S = o. (2.7) 

In these equations the unknowns are the four components of the vector 
C a , and the scalar S. The specific pressure p is a known function of / (i.e. 
of C a , by 2.4) and of S 

p = p(C a C a ,S) (2.8) 

1 That is, in general relativity where time and length have the same dimension and 
mass-energy the dimension of a length, f and therefore x have dimension of the inverse of 
a length. 
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determined by inverting the relation 2.2. 

There seems to be more equations than unknowns, but the equations 2.5. 
are not independent, because they satisfy the following identity: 

C^iVaCp - VpC a ) EE 0. 

The equation 2.7 says that S is constant along the flow lines, hence constant 
in spacetime if constant initially. We suppose, to simplify what follows be- 
low, that S is constant. Removing this hypothesis introduces no essential 
difficulty, but one has to add as new unknowns the derivatives diS, and the 
equations obtained from 2.7 by taking a derivative Vj. Our results are valid 
as they stand for a barotropic fluid, because then fj, = n(p). 



3 Bianchi equations. 

We now work in a Cauchy adapted frame, that is, a frame with its timelike 
axis orthogonal to the space slices. The spacetime metric takes then the 
usual 3+1 form 

ds 2 = -N 2 dt 2 + g^dx 1 + (3 i dt)(dx j + (3 j dt). 

The derivatives d a are the Pfaff derivatives in the coframe 6° = dt, 9 l = 
dx % + f3 l dt, that is 

d ■ d d 
dt dx^ 1 dx i 

We have written in CB.Yo 1997 the Bianchi equations satisfied by the 
Riemann tensor as a FOSH system 

Vo-R/ifc,oj + Vfc-RoMj ~ ^ hRokfij = (3.1) 

and 

^VC.u, + ^hRt.iflj — -A.il/ = VoPji - VjPoi (3.2) 

The equations 3.1. and 3.2. are for each given pair (Oj) a first order sym- 
metric system, hyperbolic relative to the space sections for the components 
Rhk,oj and Roh,oj because the matrix M° of the coefficients of the derivatives 
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8q is the unit matrix, and the matrix M i of the coefficients of the derivatives 
is indentical to M°. 

Analogous results hold for the components Rhk,ij and Roh,ij- 
To the Bianchi system is associated its Bel - Robinson energy density on 
a space slice, namely 

B = ^(|E| 2 + |H| 2 + |D| 2 + |B| 2 ) 

where E, H, D, B are the electric and magnetic fields and inductions space 
2- tensors associated with the Riemann tensor|: 

Eij = R°i t oj, Dij = -r]ihkf]jimR hk,lm , ■ 



Hij = 2^ lr lihkR hk ,oj , Bji = -N 1 r] ihk R j, hk . 
This energy satisfies the equality 

^ (|E| 2 + |H| 2 + |D| 2 + |B| 2 ) + V h {Nr) l \(E«H l:j - B i W lj } = Q(E, H, D, B) + S 

where Q is a quadratic form in E, ...B with coefficients ViV and iVK, K 
extrinsic curvature of the spaceslices. The source term S, zero in vacuum, is 

S = - ~NJ lmi r) h lm B ih . (3.3) 

The tensor p a p for a perfect fluid is given in terms of C a by 

p a /3 = {n + p)f~ 2 C a C p + -g a p{ii - p) (3.4) 

The sources of the Bianchi equations are therefore of the first order in C a , 
linear in the derivatives V a C^. 

2 Vijk the volume form, ||and V, the norm and covariant derivative, are defined by the 
space metric gy.. 
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4 Equations for VC. 



The dynamical acceleration VC satisfies the following equations obtained by 
covariant differentiation of 2.5. and 2.6, and use of the Ricci identities: 

M 7/3 = C a {V a C lP - V C ja ) + a 7/3 = (4.1) 

and 

g aP V a C^ + (/4 - l)^-V a C 7/3 + 6 7 = (4.2) 
where we have set 

C 7/3 = V 7 C /3 , (4.3) 

a 7/3 EE C y a (C af 3 — C/3 a ) + C a C\R^p X (4.4) 

by = ~R,xC x + V 7 {(/4 - l)-^-}C aP (4.5) 

The last term in 6 7 is a quadratic form in C a p whose coefficients are functions 
of the C a and S. These functions can be computed by using the identity (we 
now use S = constant) 

V 7 /i p = // p 2d 7 p, // p2 (p, S) = — . 

with by the definition of / and the identity 2.4 it holds that 

d lV = (/i + P )r 1 d,f = -(/i + v ){c x c x y i c a c ia . 

We complete the computation by using 2.1 and 2.8. 

The equations 4.1. are not independent, because they satisfy the identities 

C p M lP ee 0. (4.6) 

The equations 4.1 and 4.2. are not a well posed system. Instead of the 4x4 
equations 4.1 we consider^ the 4x3 ones: 

M 7 . = M 7i - ^M 70 = (4.7) 



3 An analogous procedure is used for the symmetrization of the Euler equations in K.O. 
Friedrichs 1969 and in Rendall 1992. 
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The terms in derivatives of C 7 a in these equations can be written in the 
following form: 

C a d a {C, % - ^C 70 ) - (d t - ^d )(C a C, a ) (4.8) 



Lemma 4.1 The system 4-2,4-7. is equivalent to a FOS (First Order Sym- 
metric) system for C ia with coefficients being functions of the Riemann 
tensor, the connection and the dynamical velocity C\, but not of their deriva- 
tives. 



Proof. The system is quasi diagonal by blocks, each block corresponding 
to a given value of the index 7. We will write the principal operator of a block 
by omitting this index. We set 

C 

Ui = C r i — — C 7 o, Uq = C a C ia (4-9) 
^0 

and we define the differential operators d a as follows: 

d = C a d a , d i = d i -^-d (4.10) 

The principal terms (derivatives of C ia ) in the equations 4.7 with index 7 
are 

doUi-diUo. (4.11) 

We have by inverting 4.9: 

C (U - C l Ui) 



C70 = 



c A c A 



The principal terms of 4.2. read, using the above formulae 
/j,'C a d a U , ., C'C' x n C^C* n 
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We introduce the positive definite (if C is timelike) quadratic form 

C { C j 



Then we find that 



9*% 



C ^ 



(4.13) 



Co C^C X 
The principal terms 4.12 are therefore 



C A C\ 



+ 9 %3 d i U j 



(4.14) 



The matrix of the coefficients of the derivatives d a in the equations de- 
duced from the system 4.2, 4.7 is 



/ - 




-d 1 


-d 2 


-d 3 \ 




-&! 


k 










-~d 2 





do 





V 


-h 








do ) 



it is symmetrized by taking the product with the 4x4 matrix 

1 

gij 

4.1 Hyperbolicity. 

A FOS system is hyperbolic with respect to the space slices x° =constant if 
the system is still symmetric when written with the usual partial derivatives, 
and is such that the corresponding matrix M* of the coefficients of derivatives 

is positive definite. It admits then an energy inequality for a positive 
definite energy relative to the space slices. This fact is of physical interest 
and also for the proof of existence of solutions of the Cauchy problem, at 
least of solutions local in time, 

In the case we are considering the matrix M° is diagonal, with positive 
elements if /i' p > and C a is timelike. The principal matrix written with 
natural coordinates is also symmetrizable. The corresponding matrix M* is 
not diagonal, and it is not obvious that it is positive definite. In fact it will 
be so only if fj,' p > 1. We will prove the following lemma. 
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Lemma 4.2 The system 4-2, is FOSE if (// > 1 and C is timelike. 

Proof. It is simpler to compute directly the energy inequality for the 
considered system: its positivity is equivalent to the positivity of the matrix 
M*. Multiplying 4.2. and 4.7 respectively by U and g lj Uj gives equations of 
the form 

^-?W,=W (4.15) 



where the F a contain only non differentiated terms. We add these two equa- 
tions, replace the operators d by the operators d and carry out some ma- 
nipulations using the expression for g % i and the Leibniz rule. We obtain 
that 

The function Qp is a quadratic form in C ja and the Riemann tensor, while 
W and T are given by the following expressions; 



T = f- 2 [(fi' p - l)U 2 + (U - C*U t ) 2 ] + g i W i U j 

The energy of the dynamical acceleration VC relative to the space slices is 
the quadratic form T . It is positive definite if fx' > 1 and C a is timelike. ■ 

Remark 4.3 The system is hyperbolic in the sense of Leray if fj,' > 0, but 
the submanifolds x° = constant are 'spacelike ' with respect to the fluid wave 
cone only if the fluid sound speed is less than the speed of light, i.e., ji' p > 1. 



5 Bel-Robinson type energy of the system. 

The Bianchi system 3.1,3.2 and analogous equations obtained by replacing 
(Oj) by (ij), together with the system 4.2, 4.7. satisfied by the dynamical 
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acceleration constitute a FOSH system if /i' p > 1. Its Bel - Robinson type 
energy (superenergy) density on a space slice is the sum of the Bel - Robin- 
son energy density of the gravitational field, and the energy density of the 
dynamical acceleration: 

£ = B + T. 

Using the expression of <9o and the mean extrinsic curvature r = g l] Kij of the 
space slices S t , whose volume element we denote by [ig t , we obtain an integral 
equality whose right hand side couples gravitational and fluid superenergies: 

f Sfi s = f Sfi s + f j {-NTS + Q G + S + Q F }i2gd6 

J St J St J to J Sg 

The scalars Qg,<S and Qp can be estimated in terms of £ and thus lead to 
a linear inequality for £, permitting the estimate of its growth with time. 
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